The CMB temperature bispectrum induced by cosmic strings 
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The Cosmic Microwave Background (CMB) bispectrum of the temperature anisotropies induced 
by a network of cosmic strings is derived for small angular scales, under the assumption that the 
principal cause of temperature fluctuations is the Gott-Kaiser-Stebbins (GKS) effect. We provide 
analytical expressions for all isosceles triangle configurations in Fourier space. Their overall am- 
plitude is amplified as the inverse cube of the angle and diverges for flat triangles. The isosceles 
configurations generically lead to a negative bispectrum with a power law decay £~ 6 for large multi- 
pole I. However, collapsed triangles are found to be associated with a positive bispectrum whereas 
the squeezed triangles still exhibit negative values. We then compare our analytical estimates to a 
direct computation of the bispectrum from a set of 300 statistically independent temperature maps 
obtained from Nambu-Goto cosmic string simulations in a Friedmann-Lemaitre-Robertson- Walker 
(FLRW) universe. We find good agreement for the overall amplitude, the power law behaviour 
and angle dependency of the various triangle configurations. At £ ~ 500 the cosmic string GKS 
effect contributes approximately the same equilateral CMB bispectrum amplitude as an inflationary 
model with |/nlI — 10 3 , if the strings contribute about 10% of the temperature power spectrum at 
£ — 10. Current bounds on /nl are not derived using cosmic string bispectrum templates, and so 
our /nl estimate cannot be used to derive bounds on strings. However it does suggest that string 
bispectrum templates should be included in the search of CMB non-Gaussianities. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 



Cosmic strings are line-like objects formed in the early 
universe P, [3, Q • They could be solitons in field theories 
with spontaneously broken symmetries [4] , singular solu- 
tions with cylindrical symmetry in supergravity theories 
}5| or fundamental objects in string theory [||. They may 
form in thermal phase transitions [H , at the end of hy- 
brid inflation 0, H, H, EH > or by tachyon condensation at 
the end of brane inflation when brane and anti-brane an- 
nihilate [ll|, [12] • If cosmic strings are added to the stan- 
dard power-law ACDM model, the Cosmic Microwave 
Background (CMB) data is fitted even better [H, Q3 
if the fraction of the temperature power spectrum (at 
I = 10) due to strings fio is about 0.1. There is there- 
fore strong motivation to develop further tests for strings 
in future CMB data, which provide the cleanest and 
best understood cosmological string 

signalsJTl]ll[i3- 
Calculations of the polarisation B-mode 

El 

Il9|| show 

that a promising line of attack for the near future is 
to look for a signal peaked between £ = 600 to 1000. 
Simulations of Planck data show that it will be sensi- 
tive down to fio ~ 0.01, and that there is no danger 
of confusing strings with inflationary tensor perturba- 
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tions [l8(. Perhaps the most characteristic signal in the 
CMB comes from the Gott-Kaiser-Stebbins (GKS) ef- 
fect HH , which is due to the gravitational lensing of 
photons passing near a moving string. This produces a 
discontinuity in the apparent temperature approximately 
proportional to the transverse velocity of the string v and 
the string tension U: 



ST ~ 8ir{GU)vT c 



(1) 



where G is Newton's constant. Given that strings move 
with a mean square (RMS) velocity of between 0.25 and 
0.36 [22], [HI we would expect to see discontinuities 6T 
of up to about 1 fiK. Recent calculations of the inte- 
grated Sachs- Wolf CMB signal at small angular scales 
have shown that the angular power spectrum it produces 
decreases slowly, approximately £~°- 9 at high multipole 
moment I [24j]. If /io — 0.1 strings should dominate 
the adiabatic temperature power spectrum for I > 3000, 
and remain above the thermal Sunyaev-ZePdovich effect. 
This fraction corresponds to GU ~ 0.7 x 10 -6 for strings 
in the Abelian Higgs model [Hj]. The slow power-law 
decrease in the power spectrum is quite close to a predic- 
tion of the small scale power spectrum using the string 
correlation functions in a Gaussian approximation [25| . 
which gave a temperature anisotropy power going as £ . 
The non-Gaussian nature of the maps exhibited in [24j |. 
as for example exhibited by the skewness of the distri- 
bution of the 1-point function, immediately motivates an 
attempt to generalise the calculation to higher-order cor- 
relators, and in particular the 3-point function or bispec- 



2 



trum, which will become increasingly well characterised 
by future CMB data [H EE El El- This paper reports 
on the results of this calculation. It is found that the 
bispectrum for an isosceles arrangement of wavevectors 
k in Fourier space is generically negative and decreases 
as |fc|~ 6 - It is proportional to (GU) 3 and hence to the 
power spectrum raised to the 3/2 power. The scale of the 
non-linearity parameter /nl, defined by dividing the bis- 
pectrum by the square of the power spectrum, is therefore 
potentially large, going as (GU)^ 1 . 



where we have defined 

e = 8ttGU, (8) 

and A = 1,2 with implicit summation on repeated in- 
dices. It is now clear that the power spectrum, bispec- 
trum, and higher order correlators can be evaluated in 
terms of correlation functions of the string network, as 
projected onto our backward light cone. In the next sec- 
tion we will introduce the relevant correlation functions 
and discuss their important features. 



II. ANALYTICAL BISPECTRUM 



B. String correlation functions 



A. Gott-Kaiser-Stebbins effect 

We work in the flat sky approximation, and define 
transverse coordinates x measured in radians. The wave 
number k is related to the multipole moment £ by [29|, [30j 



m + i) 



(2) 



The anisotropy power £(£ + l)C? is then approximately 
equal to k 2 \5Tk\ 2 , where STf. is the Fourier transform of 
the temperature fluctuation, 



dx8Te ikx . 



(3) 



We will also define Q(x) — 5T(x)/T CM b- String space- 
time coordinates will be denoted X^(t, ct), where r and 
a are timclike and spacelike worldsheet coordinates re- 
spectively. In the temporal gauge the X a = r (where the 
worldsheet time is identified with the back gro und time 
coordinate) the (corrected) GKS formula is [25| 



V 2 6 = SttGU / da 



x _(JiJlx ■V6^(x-X), 
(X-p) 

(4) 

where p M = p^/E, and worldsheet variables are evaluated 
at the retarded time t r = t + z — X 3 (a, t t ), when the CMB 
photons (taken to be moving in the — z direction) pass 
the string. The expression is greatly simplified in the 
light cone gauge, 



(5) 



where X^ — X°±X 3 . The time parameter r then labels 
the intersections with a set of null hyperplanes with the 
worldsheet: the photon geodesies = + Ap M form 
just such a set. Then we find 

V 2 9 = -8itGU J daX ■ V5 {2 \x - X), (6) 

where worldsheet quantities are now evaluated at r = 
x + = t + z. In Fourier space the equation becomes 

-k 2 e k = iek A J daX A (a)e lkX ( a \ (7) 



We denote the transverse coordinates of the string, 
X A (a). The basic two point functions are 



(X A (a)X B {a')), (X A (*)X B (*')), (X A (a)X B (a')), 

(9) 

where the angle brackets denote an average over an en- 
semble of strings. The starting assumption is that the 
string ensemble is well-approximated by a Gaussian pro- 
cess: that is, all correlators can be calculated in terms of 
the two point functions. We now make some assumptions 
about the ensemble: (i) rotation, reflection and trans- 
lation invariance of the transverse coordinates; and (ii) 
worldsheet reflection and translation invariance. Then 
there are three independent correlation functions: 



(X A (a)X B (a')) 
(X A (a)X B (a')) 
(X A (a)X B (a')) 
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6 AB V(a-a% 
S AB T(a - a'), 
8 AB M 1 (a-<r'). 



A fourth, 



M 2 (a) = (X A (a)X B (0))t 



(10) 

(11) 
(12) 

(13) 



vanishes because of the symmetry X 1 «-> X 2 (this point 
was overlooked in Ref. [25|]). The functions V and T are 
symmetric in their argument, while M\ is antisymmetric. 
The forms of the correlators are sketched in Fig. Q] (see 
also [111]). For later convenience two other correlators 
will be defined: 



[X(a)-X(a')] 

per p(T 

= d<7 1 da 2 T(a 1 - cr 2 ), 

J a' J a' 

n(a-*') = ([X(a)-X(a'))].X(*') 
= f daiMi(ax - a'), 

J a' 



(14) 
(15) 
(16) 
(17) 



There are integral constraints arising from the fact that 
the average velocity of the strings vanishes on large scale: 



daV(a)^0, daM^a) 



0. 



(18) 
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FIG. 1: Sketches of string correlation functions V(a) (velocity- velocity) , T(a) (tangent-tangent) and iWi(er) (velocity-tangent), 
as a function of the string worldsheet spacelike separation a, defined in Eqs. (|10[) to (|12[) 



The important asymptotic forms are: 





r v 2 a^o 

a — ► oo ' 


(19) 


T(a) - « 


f Pa 2 a -» 
£<7 cr — > oo ' 


(20) 


n(<r) -> < 


a — > oo 


(21) 


where we have defined four parameters 




i = 


r'(oo), 


(22) 


v 2 = 




(23) 


p = 




(24) 


CO = 




(25) 



The equations of motion in a FLRW background follow 
from the Nambu-Goto action 



S = -U J drdcV-7 



= -U J drda^(X u )V -X 2 X 2 + (X ■ X) 2 , (28) 

where 7 is the determinant of the induced metric along 
the string worldsheet, and a(X°) is the scale factor. With 
the standard gauge choice X ■ X = 0, we find 



\e a/ e Oa \e Oa 

1 da 1 dX° 1 dJf ^ 1 da • 2 _ 

" adX^l~da~7~da~ + °adX^ ~ °' 



(29) 



where e = y — X 2 /X 2 . The light-cone gauge choice 
X + = r produces as the equation of motion for X + 



The correlation length £ is the projected correlation 
length on the backward lightcone, P is the mean square 
projected tangent vector (of order unity), v 2 is the mean 
square projected velocity (again of order unity), and cq 
is the correlation between the projected velocity and cur- 
vature. 



C. Light cone gauge equations 



2H(X° + X 2 ) = 0, 



(30) 



where TL = dlna/dX is the conformal Hubble parame- 
ter. The equation for the transverse components is 



1 / - 2\ • 1 d fldX 
X + 2TL-[X )x--—[~ — 

e A V / £ oa \ £ oa 



2H-t(X ■ X)X = 0. 



(31) 



In Minkowski space, the Nambu-Goto action leads to 
the equations of motion and constraints 



X t *-X>* = 0, X 2 +X 2 = 0, X-X = 0. (26) 

The light-cone gauge consists of choosing r = X + = 
X° + X 3 . Hence 



X~ = X A X A + X A X A , X~ = 2X A X A . (27) 



In a FLRW background, assuming that (X ) is constant, 
and neglecting higher-order correlations between TL, X 
and X, we find 



OS 



(32) 
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where we have defined ds — edo, and where TL is the 
conformal Hubble parameter averaged over the string en- 
semble at fixed t — X + , which will select the value where 
there is most string, i.e. at decoupling. 

If one we assume that the string ensemble is approxi- 
mately Gaussian in X and X/e, it is not hard to show 
that the right hand side reduces to 



d 2 x 

ds 2 



X 



H\ (X 



X X 



(33) 

The last term vanishes, so the cross correlator simplifies 
to Hv 2 P which is positive. From Eq. (|25p . one gets that 
in a FLRW background, a string network should exhibit 
Co > 0. We see also that it vanishes in Minkowski space, 
which can be viewed as a consequence of time reversal 
invariance. As shown in the next sections, the temper- 
ature bispectrum vanishes when the correlator II does, 
and so the generation of a bispectrum by strings simply 
requires the breaking of time reversal invariance, as it is 
in a FLRW background 1 . 



D. Temperature power spectrum at small angular 
scales 



The power spectrum at small angular scales, where 
the GKS effect [20l l3lT | is held to be dominant, was cal- 
culated in Ref . |25| , and the calculation is recapped here 
for completeness. The Fourier transform of the temper- 
ature fluctuation given in Eq. ^ reads 



(34) 



and hence is given by Eq. ([7]). The power spectrum is 
defined by 



(e fcl e fc2 > -p(fc 1 )(2^) 2 <5(fe 1 + A ;2 ) 



(35) 



With our conventions we need a box of formal area A = 
(2ir) 2 5(0) to express the power spectrum 



2 k A k B 



xe 



ik-[X(a)-X(a')} 



(36) 



With our assumptions about the string correlation func- 
tions, the ensemble average can be reduced to 



i 



dodo' 



-eiPio-o 1 ) 

2 v ' 



V(o - o') 



-fc 2 r(tT-<r')/4 



(37) 



We now derive the asymptotic behaviour of P(k) as (fc£) 
gets large. The contribution to the power spectrum from 
the mixed correlator M\ can be shown to be subdominant 
at high k [25| , and so we need examine only the first term 
in the power spectrum of Eq. ([37]) , 



P(k) 



! ^ 2 Jdo + da.V^)e 



fc 2 n>_)/4 



where o± = o ± o' . For fc£ 3> 1, we find 



k 2 P(k) ~ e 2 0F 



v 2 1 



(38) 



(39) 



where L is the total transverse light-cone gauge length of 
string in the box of area A. 

The power spectrum given in Ref. [12] was consis- 
tent with k , and the amplitude was surprisingly close 
for such a crude estimate. Fraisse et al. [24| have a 
slightly different small-scale angular power spectrum: 
k 2 P(k) ~ k~ p with p ~ —0.89. This can be explained 
if r(cr) ~ o 2 l p on the relevant scales. This correlation 
function controls how far on average one moves in the 
transverse coordinates as one moves along the string: 
p = 1 would correspond to straight lines, while p — 2 to 
a Brownian random walk. A power less than one is sug- 
gestive of a cloud of zero-dimensional objects along the 
string worldsheet which may be the signature of small 
loop production. 



E. Temperature bispectrum from strings at small 
angular scales 

In the flat sky approximation the three point temper- 
ature correlation function or bispectrum is defined as 

(e fcl e fc2 e fc3 ) -s(fc 1 ,fc 2 ,fc 3 )(27r) 2 ( 5(fc 1 +fc 2 +fc3). (40) 

We again need to normalise by a formal area factor A = 
(2tt) 2 <5(0) to obtain an expression in terms of a string 
expectation value 



B(k 1 ,k 2l fc 3 ) = ie 3 -rS AA S BS S C c ijAr§ 

x J do.dozdo-^XfxixCe^^x*), 

(41) 

with = X A (o a ), a, be {1, 2, 3}, and k x + k 2 + fc 3 = 
0. With the Gaussian assumption, the ensemble average 
of the string observables is lengthy but straightforward. 
Defining 



C 



ABC — vA yB V C 



XtX*X^ , D = 5 ab k a ■ X b , (42) 



1 We thank B. Wandclt for explaining this point to us. 



we have 



(C ABC e lD ) =i(C ABC D)e-^ D2 \ (43) 
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and hence 
B(k 1: k 2 ,k 3 ) = -e 



we have 



3 1 h A k 2B k 3c 

A l,2jL2L2 
\ *"2 3 



/ da 1 da 2 da 3 (C ABC Z?) e -^ D2 >. 



(44) 



We first evaluate the correlator 
(C ABC D) = lx?X*X°k a ■ X a 

= (x A x B ){x%k a ■ x a ) + (x^x A )(x£k a ■ x a ) 

+ (X B X^)(X A k a -X a ). 

The velocity correlators have already been given in 
Eq. ([TO]) . The mixed correlators are most easily eval- 
uated using variations on the theme 



k a -X a = k,- (X, - X 3 ) + k 2 ■ (X 2 - X 3 ) 
= fei ■ -^13 + k 2 ■ X 23 , 

where X a b = X a — X^. For example, 

[x A k a ■ x a ) = k? (x A xg) + fcf (x£x& 



(46) 



(47) 



Using the definition of the correlator H(a) in Eq. (fT7|) , 



x A k a -x a ) = i^n(a 21 ) + ^ 3 4 n(<T 31 ), (48) 

where a a (, = cr a — crj. Substituting for the velocity corre- 
lators 



(c abc d) = \s AB [fc 1 c n( ( r 13 ) + fc 2 c n( ( T 23 )] v(* 12 ) 

+ ±5 CA [kf n(oia) + fc 3 B n( ( T32)] V(a 31 ) 
+ -8 BC [k A U(a 21 ) + kfn(a 31 )] V(a 23 ). 



Now we evaluate 

i_> ■ / . „ -r x2 



(49) 



(D 2 ) - ((fc a ■ X a ) 2 ) = (fa ■ X 13 + k 2 ■ X 23 

= ^?r(a 13 ) + \k 2 2 T(a 23 ) + fc x ■ k 2 (X 13 ■ X 2 



(50) 



It can be shown that 
1 



X 13 ■ X 23 ) = - [T(a l3 ) + T(a 23 ) - T(a 12 )] , (51) 
hence 

(D 2 ) = -l[fc x • k 3 T(a 13 ) + k 2 ■ k 3 T(a 23 ) 

+ k 1 -k 2 T{a l2 )]. 

Note that this expression is symmetric under the inter- 
change of any pair {cajO^} and {k a ,kt,}. Contracting 
Eq. |gl|) with k A k 2 k 3 , and defining 



., 1 1 

' d(Tidcr 2 d(T3 



B(^k 2 ,k 3 ) = -e- w[m 
X |k 12 [Ki 3 IT(cri3) + K 2 3ll(cr 23 )] V(o-ia) 

+Ki 3 [Ki 2 n(cri 2 ) + K 23 n(cr 32 )] 7(0-31) 

+K 23 [Ki 2 II(cr 2 i) + Ki 3 II((T3i)] V(cr 23 ) 



x exp |-- [Ki 3 r(cri 3 ) + K 23 r(er 23 ) + Ki 2 r(ai 2 ) 

. (54) 

We perform the integrations over the string coordinates 
a a in the Appendix 1X1 where it is found 



B{k 1 ,k 2 ,k 3 ) = -e 3 irc 



^ 1 



1 



3/2 



Kab — k a ■ k\) 



(53) 



t 4 A ^ k\k 2 k 2 

kfn 23 + fc|«3i + k$Ki2 
(K 23 K 31 + Kl 2 K 3 l + Kl 2 K 23 ) 

"(55) 

This is our primary expression for the bispectrum in- 
duced by the GKS effect in cosmic strings. It is propor- 
tional to (Gt/) 3 , and therefore goes as the 3/2 power of 
the power spectrum. The factor L£/A is a geometrical 
factor of order unity, as the projected string length per 
unit area is of order the projected correlation length £ 
(unless there are a large number string networks which 
do not interact with each other). The factor in curly 
brackets is geometrical in Fourier space, depending on 
the relative lengths and angles of the k a . The overall de- 
pendence on angular quantities is therefore £ -2 fc~ 6 . The 
dependence on string correlators appears through f, w 2 , 
and Co defined in Eqs. (j2"2"|) to (f2"5)) . This last is an inter- 
esting quantity as it is not time-reversal invariant: hence 
the bispectrum could vanish only if the string network is 
time symmetric. In an expanding universe, the existence 
of an asymmetry is ensured by decay of the string net- 
work. It was argued at the end of Section III CI that cq 
should be positive. 



F. Symmetrical triangle configurations 

Current analysis of the CMB temperature bispectrum 
focus essentially on a particular local model of primordial 
non-Gaussianity as well as specific configurations of the 
primordial bispectrum wavenumbers. From Eq. (|55p. we 
can easily estimate the string induced bispectrum for var- 
ious symmetrical triangle configurations in Fourier space. 
Let us stress that Eq. (f55|) is the CMB temperature bis- 
pectrum as opposed to a primordial bispectrum which 
is the three point-functions of the primordial Newtonian 
potential, usually of inflationary origin. The latter has 
still to be evolved through the CMB transfer functions. 
As a result, even for identical wavenumber configurations 



6 



40 
30 
S 20 

S 10 
m 

| -10 
*t 

+ -20 
-30 
-41) 





0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 U 
6 (rad) 



FIG. 2: Angular dependency of the isosceles bispectrum as 
a function of the angle 9 in between the wavevectors fei and 
&2. The particular values # cq = 7r/3 corresponding to the 
equilateral configuration as 9q at which the bispectrum van- 
ishes are represented. Notice the divergences for flat triangle 
configurations at 9 — > (squeezed) and 9 — > ir (collapsed). 



2. Squeezed triangles 

This is the case of one of the sides of the triangle van- 
ishing &3 — > kO, with the opposite angle 9 — > 0. The 
angular factor in Eq. (|58p simplifies to 1/9 3 and 



/,■/,■» 



limAsoCM) 



9^0 



-6 7TC 



v 2 L£ 1 1 
FT4 |2^i^3 



(60) 



The bispectrum is therefore negative and appears to di- 
verge in the squeezed limit. In practice, it is not possible 
to observe such a divergence as the maximum value of the 
wave vector is bounded by the detector resolution, and 
the minimum by the map size. Let us also recap that the 
above calculation breaks down at low wave numbers, and 
thus for the too small values of 9, since the GKS effect 
would no longer be the dominant source of temperature 
anisotropies. It is unclear at what £ this happens, but 
we note that there is no sign of an l~ x power law in [33| . 
so the GKS effect may be a subdominant contribution to 
the cosmic string power spectrum for £ < 1500. 



in Fourier space, the resulting temperature bispectra are 
not the same. 



3. Collapsed triangles 



1. Isosceles triangle 

We now consider isosceles triangle configurations in 
Fourier space such that 



Ifcxl 



\ko 



k, Ifea 



2k sin 



2' 



(56) 



where 9 denotes the angle between the wavevectors ki 
and k,2- The cross scalar products simplify to 

Ki2 = k 2 cos 6, k 23 = k 3 i = 2fc 2 sin 2 (6»/2), (57) 

and the isosceles bispectrum reads 

v 2 LS, 1 1 + 4cos6»sin 2 (6»/2) 



Biso(k, 



-e 3 7rc 



i A A £2 fc 6 



(58) 



Notice that for 9 = tt/3, we obtain the peculiar case 
of an equilateral triangle [see Eq. ([64]) ]. In Fig. [2] we 
have plotted the angle dependency of the isosceles bis- 
pectrum. Such a configuration diverges in the two flat 
triangle limits for which either 9 — > or 9 — > n. Both 
of these configurations are therefore better suited than 
the equilateral one to characterize the strings and are 
discussed in the next sections. Notice also the change of 
sign which occurs for the angle 



9n — 2 arccos 



\/3- \/3 



(59) 



For 9 < 9o, one has -B; so < (as for example in the equi- 
lateral configuration) , whereas for 9 > 9 a the bispectrum 
S iso > 0. 



There is however another flat triangle configuration for 
which all wavenumbers scale similarly. For 9 — > n, we 
have a collapsed triangle with k\ — ki = k and k% — > 2k. 
Denoting by 



if = 7T — (/, 
the isosceles bispectrum reduces to 



(61) 



B, 



lim Bi S o(M) 



-e 37tcq 



v^LiJ- 1_ 



( 62 ) 

Notice the divergence but with a positive bispectrum. 
Again, such a divergence is not observable in a realistic 
situation due to the finite detector resolution. Neverthe- 
less, flat triangles should provide the best framework for 
the search of a string bispectrum signature. In fact, the 
change of sign in between the squeezed and collapsed tri- 
angles is also of interest to improve the signal to noise 
ratio. Substracting these two configurations with the ap- 
propriate angles should enhance the string signal over the 
noise. 



4- Equilateral triangles 

For completeness, we give the bispectrum in the par- 
ticular case 9 — 7r/3, where the k a are arranged in an 
equilateral triangle: 



Ifcil 



lfe 2 



\ko 



(63) 
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One has K a b = k 2 /2, and the factor in the squared brack- 
ets in Eq. ([55)) becomes 4/V3. As the result, 



Bfefefc = B cq (k 1 ,k 2 ,k 3 



Air v 2 L£ 1 



(64) 



With co > 0, such a configuration produces a nega- 
tive bispectrum decaying as 1/k 6 but with an overall 
amplitude significantly smaller than the collapsed and 
squeezed triangle configurations. 



III. NUMERICAL ESTIMATE 

In this section, we directly compute the three point 
function in Fourier space from a set of 300 CMB temper- 
ature maps induced by cosmic strings in the flat sky ap- 
proximation. These maps have been obtained along the 
lines of Ref. [24 ] from Nambu-Goto numerical simula- 
tions in FLRW space-time. As detailed in this reference, 
they are obtained from the GKS effect using Eq. (@|, 
and are therefore valid on small angular scales only. In 
the next sections, we briefly recall the numerical method 
used to generate the maps and then discuss our bispec- 
trum estimator. 



A. Nambu-Goto numerical simulations 



procedure introduces one additional parameter which is 
the redshift at which we start the simulations z\. It has 
been set to the last scattering surface, namely Z\ = 1089, 
and in a flat "Lambda-Cold-Dark-Matter" (ACDM) uni- 
verse, using fiducial values for the density parameters 
compatible with the three-year Wilkinson Microwave 
Anisotropy Probe (WMAP) data (37j, this corresponds 
to a numerical comoving box of L s im ~ 1.7 Gpc. Such 
a size subtends a angle of 6f ov ~ 7.2° in the sky. The 
Nambu-Goto simulations and the associated CMB maps 
therefore depend on only two parameters: the string en- 
ergy per unit length U, and the initial correlation length 
£ c . The dependence on £ c should drop out at late times, 
as the network is believed to approach a self-similar or 
scaling configuration. However, in a real Nambu-Goto 
simulation, there is still some non-scaling structure evi- 
dent on the smallest length scales [24] . As an illustration 
example, a typical CMB temperature anisotropy map is 
represented in Fig. [31 together with the seeding strings 
projected on our past light cone. We have used our nu- 
merical simulations to create 300 statistically indepen- 
dent temperature maps from which one can construct a 
bispectrum estimator. 



B. Reduced bispectrum estimator 

1. Scale convolution method 



Our FLRW numerical simulations are based on an im- 
proved version of the Bennett and Bouchet Nambu-Goto 
cosmic string code [24| H3, H|[ • The runs are performed in 
a comoving box with periodic boundary conditions and 
whose volume has been scaled to unity. The horizon size 
e2h is a free parameter controlling the initial string en- 
ergy within a horizon volume. For the simulations we 
performed, dh a ~ 0.185. The string network is assumed 
to come from Vachaspati-Vilenkin initial conditions for 
which the long strings path is a random walk of corre- 
lation length £ c , plus a random transverse velocity com- 
ponent of root mean squared amplitude 0.1 [36l ]. These 
parameters are set as in Ref. [24| to minimize the re- 
laxation time of the Vachaspati-Vilenkin string network 
toward its stable cosmological configuration. The tem- 
perature maps are then produced according the GKS ef- 
fect generated by the strings intercepting our past light 
cone, using the Eq. (I) 
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iGU 



X n: 



(u ■ k) e 



-ik-X 



e- T ed<T, (65) 



where t(X) is the optical depth to the position of the 
string, along the line of sight n, and 



u = X 



(ft ■ X 1 ) ■ X' 
1 + n-X 



(66) 



The cosmic string simulations are used to compute u, 
based on the string trajectories X. The map generation 



The bispectrum computations use the scale convolu- 
tion method introduced in Ref. 1 26] and applied to the 
flat sky approximation in Ref. [38|. This method relies 
on the choice of unity window functions in Fourier space 
W u {k) peaked around a particular wavenumber u. Defin- 
ing 



e„(a:) 



dfc 



Q k W u (k)e 



-ik-x 



(67) 



one can construct an estimator of the three point function 
in Fourier space by remarking that 



6 fcl (x)Q k2 (x)Q ka (x)dx 



(27T) f 



(68) 



x W kl (p) W k2 (q) W k3 (k) (2tt) 2 S( P + q + k). 



For peaked enough window functions, & k remains con- 
stant over the window functions width and we construct 
our reduced bispectrum estimator as 



1 



kxk 2 k 3 



s 



e kl (x)Q k2 (x)Q k3 {x)dx) . (69) 



where the function is the flat sky equivalent of the 
inverse Wigner-3j symbols and reads 



k±k 2 k 3 



^Wk 1 (p)W k2 ( q )Wk 3 (\p + Q\)- (70) 
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FIG. 3: Typical CMB temperature map on a 7.2° field (res- 
olution of 9 ICB = 0.42', np ix = 1024 2 ) obtained from Nambu- 
Goto cosmic string simulations (upper panel). The lower 
panel traces back the strings projected on our past light cone. 



As noted in Ref. [38|, S is of geometrical nature and 
needs to be computed only once. However, at small an- 
gular scales, S is generically a four-dimensional integral 
whose computation can be time consuming for the large 
wavenumbers. For thin enough window functions, it is 
nevertheless possible to derive an analytical approxima- 




1000 

i 



FIG. 4: Exact versus approximated geometrical S^t factors 
for a map with 1024 2 pixels and various choices of the window 
function with w (w = 2n corresponds to one Fourier mode per 
bandwidth) . For large values of w, the deviations for the small 
multipoles come from the averaging effect. For small angles 
we are interested in, Eq. (|72p provides a good approximation 
to the exact expression (|70ll . 



tion. Assuming that W u (k) = 1 for u — w/2 < k < 
u + w/2, and zero otherwise, for small enough width w 
one has 



W u (k) ~ w5(k-u). 



(71) 



With fci, k-2 and ks forming a triangle, Eq. ([70]) can be 
exactly integrated and one finds 



Ak 1 k 2 k 3 



(h + k 2 y 



fc 3 



fc 3 



(fci - k 2 y 



(72) 

As can be checked in Fig. 2J the analytical approxima- 
tion of Sfc™fc 2 j. 3 given in Eq. (|T2")) is particularly accurate 
for the large wavenumbers. This approach ends up be- 
ing numerically convenient since it requires only three 
Fourier transforms to compute the 9fc(cc), together with 
an integration over all map pixels. 



2. Statistical averaging and numerical tests 

From the previous discussion, the bispectrum is ex- 
tracted and averaged over the 300 statistical indepen- 
dent maps and for the various triangle configurations dis- 
cussed in Sect. Ill F] Our results are presented in the next 
section. The variance of the bispectrum over the differ- 
ent maps is used to provide an estimate of the statistical 
errors associated with our approach. Notice however, as 
discussed at length in Ref. [2J], the Nambu-Goto string 
simulation do also lead to some systematic errors due to 
the non-scaling structure. This can be greatly reduced 
by eliminating some of the small loops from the network. 
In order to check our numerical bispectrum estimator, 
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FIG. 5: The reduced bispectrum computed by the scale convo- 
lution method and averaged over a hundred of non-Gaussian 
synthetic ACDM temperature maps whose probability distri- 
bution function is given by Eq. (|73[1 . The red dashed line 
is the expected non-Gaussian signal. The missing power for 
low multipoles comes from the cutoff at £ w ~ w/8t ov ~ 1050 
introduced by the choice of large w values for the window 
functions. 



timator 2 at the price of losing information for the lower 
multipoles and those close to the Nyquist frequency. The 
second test we have performed is to integrate the cosmic 
string bispectrum over all possible triangle configurations 
such as to check that we indeed recover the skewness of 
the string temperature anisotropies. Averaged over the 
300 string maps, we find the mean sample skewness to 
be negative 



9x 



_ I (e-e) 3 



-0.22 ±0.12, 



(75) 



where the brackets stand for the mean over different re- 
alisations while the bar denotes averaging on each map. 
The variance itself averages to 



(e-e) 2 ) ~ (150.7 ± is) (Guy 



(76) 



The quoted errors are statistical and refer to the square 
root of the variance between the different realisations. 



C. Flat triangle configurations 



we have performed two tests. The first was to generate 
a set of synthetic non-Gaussian CMB temperature maps 
from the non-Gaussian probability distribution 



The previously described method has been applied to 
extract the isosceles bispectra in the flat triangle config- 
urations for which the amplitude is maximal. 



iV(e) 



3 -e 2 /(2<r 2 ) 



27TCT 



H, 



e 

V2a 



< 73 ) 

where H3 stands for the third Hermite polynomial. As 
shown in Ref. [1^, H3] , such a statistics leads to a reduced 
bispectrum whose amplitude is given 



y/ Ce 1 Ci 2 Ce 3 n p i x ^jV 2 ' 



(74) 



where /i 2 = cr 2 (l + 6a 2 ) and fi 3 = (2cr 2 ) 3 / 2 a v /3(l - a 2 ) 
are the second and third central moment, respectively. 
It is not difficult to show that /x 3 / ' [Sj 2 i s maximised for 
a 2 = (7— V43)/6. In Fig.[5J we have plotted the bispec- 
trum obtained from the scale convolution method aver- 
aged over a hundred of such non-Gaussian ACDM tem- 
perature maps together with its analytical expectation. 
At small scales, and up to the statistical errors, the non- 
Gaussian signal is recovered with the right amplitude. 
Notice that the loss of power for the low multipoles is sim- 
ply an artifact coming from the choice of a large window 
function width (w — A2n). This necessarily introduces 
a lower frequency cutoff around i w ~ w/0[ ov ~ 1050, 
and also close to the Nyquist frequency (not visible on 
the plots). As discussed in the following, increasing the 
values of w reduces the variance of the bispectrum es- 



1. Squeezed triangle 

Let us first consider the triangle configuration in the 
squeezed limit such that 



k x =k 2 = k, 



k9, 



(77) 



that we refer to as "U0" . The mean value of bug over the 
different maps have been plotted in Fig. [6] and for various 
values of the squeezed angle 9. As found analytically, 
the bispectrum is negative while the overall amplitude is 
found to be enhanced by a factor 1/9 3 . We have used a 
window function width equals to w — IOit for this plot. 
However, due to the finite field of view, we cannot go to 
very small squeezed angle without truncating the lower 
modes. From Eq. (|77|) . the lowest multipole achievable is 
therefore given by 



t ■ = 



w 



7fovC 



(78) 



where the window function cutoff has also been included. 

In order to determine the power law behaviour of the 
bispectrum , we have plotted in Fig. the mean value 
and standard deviation over the 300 different maps of 



2 Simply by increasing the number of equivalent triangle configu- 
rations in Fourier space probed by the window functions. 
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FIG. 6: Squeezed angle dependency of the mean squeezed bis- 
pectrum [£(£ + 1) / (2n)] s/2 bue- The spurious plateau (dotted 
lines) for the lower multipoles comes from the cutoff associ- 
ated with a given squeezed angle 8 together with the window 
function width w and occurs at £ m i n = w/(9t ov 8). 



FIG. 8: Rescaled mean squeezed bispectrum 9 3 [£{£ + 
l)/(27r)] 3/,2 6«e showing the I/O 3 dependency. Notice the large 
window function width w = 30n chosen to reduce the statisti- 
cal errors and pushing up the numerical multipole cutoff l m in- 



| 6 = 0.20 rad | 




Fig. the standard deviation of our estimator becomes 
an order of magnitude bigger than the bispectrum itself 
for £ > 16000. 

The power dependency of Eq. ([79)) is remarkably close 
to the analytical expectation. In fact, the non-integer 
value for the power may be expected due to the fractal 
properties of the Nambu-Goto strings [23| . As shown in 
Ref. [2J| , the power spectrum itself was found to behaves 
as 



1(1 + l)C e cx r p with p = 0. 



(80) 



From Eqs. (|T5|) and ([50")) , we the find the analytical power 
law behaviour of the ratio 



FIG. 7: Mean value and standard deviation of the squeezed 
bispectrum [£(£ + 1)/(2tt)] 3/2 b U e for 8 = 0.2 radians. The 
window function width has been increased to w = 30-7T to 
reduce the variance pushing up the lower multipole cutoff to 
imin = w/(8f ov 8) __ 3750. The dashed line is the best power 
law fit whose power has been set to the one obtained from the 
equilateral configuration [see Eq. (|84|l ]. 



the squeezed bispectrum for the particular angle 6 = 0.2 
radians. The variance has been reduced as far as we could 
by using a large window function width w = 30n. This 
explains the low multipole cutoff at a quite large value of 
£ m i n ~ 3750. A power law fit against the mean numerical 



estimator for £ 
gives 



> ^min and truncated at l„ 



16000 



[lit + 1)] 3/2 b ue oc r q with 9 = 2.82. (79) 

We do not quote error bars since this value is relatively 
sensitive to the choice of £ ma x and, as can been seen in 



oc r {q - 3p/2) with ? --p=1.49. (81) 

^3/2 e>1 2 

The dashed blue line in Fig. [7] is the best power law fit 
£~i. 

The I/O 3 dependency of the squeezed bispectrum is 
illustrated in Fig. [51 Again, to reduce the variance, the 
width has been set to to = 307r and we have plotted 
9 3 [£(£ + l)/(2ir)] 3 / 2 b ue for the four 9 values of Fig. © 
As expected, all rescaled mean values match above their 
respective multipole cutoff £ m i n [see Eq. (175)) ]. 

Concerning the overall amplitude, it can be evaluated 
around the minimum variance multipole and we find at 
£ = 5000 

[£{£ + l)/(2^)] 3/2 b m ~ (-2.7 ± 1.4) x lO" 3 (?pj " . 

(82) 

To conclude this section, the squeezed bispectrum is 
found to follow the analytical results of Sect. [TT1 up to 
the non-integer power law behaviour with respect to the 
multipole moment which may be interpreted as a fractal 
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FIG. 9: Mean collapsed bispectrum [t{t + l)/(2w)] 3/2 be v>v> . 
This configuration is extremely sensitive to the window func- 
tion properties and statistical errors dominate. The change 
of sign for small values of ip is nevertheless observable. 



effect associated with the Nambu-Goto string small scale 
structure. 



bispectrum being of smaller amplitude that the squeezed 
one, its numerical determination is quite challenging. It 
is however clear from this plot that bui < in agreement 
with the analytical results of Sect. [TTJ 




FIG. 10: Mean and standard deviation of the equilateral bis- 
pectrum [£(£ + l)/(27r)] 3/ ' 2 but over our 300 string tempera- 
ture maps and obtained for window function width w = 307r. 
The blue dotted line is the power law fit of Eq. (l84t . 



2. Collapsed triangles 



As discussed in Sect. Ill F~3l the collapsed triangles cor- 
respond to the limit ip — > for the isosceles configuration 
and to the mode values 



k. 



2k. 



(83) 



The scale convolution method is hardly applicable in this 
case since all collapsed triangles fitting in the window 
function width will be indistinguishable. Although this 
averaging effect allowed us to increase the statistics for 
the squeezed triangle configurations, for the collapsed 
ones one cannot longer assume that the bispectrum is 
constant over the window function, especially since it di- 
verges in the limit ip — > 0. As a result, w has been kept to 
a low value (w — 3tt) and the modes have been binned to 
reduce the statistical errors. As can be seen in Fig.[9l our 
numerical estimation does not allow to accurately traces 
the bispectrum down to the small values of ip. Neverthe- 
less, we do observe the change of sign for ip — ► albeit 
at a smaller values than the one predicted analytically. 



D. Equilateral configurations 

The previously described method has also been applied 
to extract the equilateral bispectrum but- In Fig. [TU1 
we have plotted its mean value over the 300 different 
maps as well as its standard deviation. To reduce the 
variance of the estimator we have chosen w = 30n which 
gives a low multipole cutoff £ w ~ 750. The equilateral 



A power law fit against the mean numerical estimator 
over the range of multipolcs £ w < £ < £ max , with £ max — 
16000 gives the same power as for the squeezed triangle 
configurations, namely 

[£(£+ l)f /2 b m oc^ r q with g = 2.82. (84) 

Concerning the overall amplitude, we find at I = 1000 
[£(£ + 1)/(2tt)] 3/2 b m =s (-0.53 ± 0.3) {GU)\ (85) 

which is again of the same order of magnitude as our 
analytical estimate. Using the power law fit, such an 
amplitude can be extrapolated up to I — 5000 and com- 
pared to Eq. (|82p. In particular, we recover, up to the 
statistical errors, the analytical result 



6 3 bm 



4 

7s 



(86) 



For completeness, the above ratio has been plotted for 
6 = 0.2 radians in Fig. El 

As can be seen in Fig.[l0l but also in Fig. [71 the bispec- 
trum decays faster than its standard deviation at small 
scales. This can be understood by using the nearly Gaus- 
sian approximation in which the non-Gaussian trispec- 
trum is neglected. Under these assumptions, the variance 
of the bispectrum estimator, for the equilateral configu- 
rations, is given by [H, [4l[ 



4 
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S. 



(w) 



c 



(87) 
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FIG. 11: Ratio of the mean bispectra but I '{d 3 bue) obtained 
for — 0.2 radians. Up to the statistical errors, the expected 
value 4/V3 is recovered for £ > t m i U - 



Measured 
Gaussian approx. 



FIG. 12: Standard deviation of the bispectrum estimator over 
the 300 string temperature maps (black solid curve) compared 
with the nearly Gaussian random expected standard deviation 
obtained from Eq. 



Using Eq. (|72|) . one gets 



2tt 



'fOY 





27 


1/4 


\i{l + l)Ci 


3/2 




U{l + \)_ 




2tt 





(88) 

As expected, the variance is reduced by a factor w 3 
coming from the different numbers of equivalent trian- 
gle configurations probed by window functions of width 
to. From Eqs. ((50)) and ([55)1 . one gets <r b oc £ -4 - 83 which 
vanishes slower than the bispectrum estimator at large 
multipoles. In Fig. we have plotted the standard de- 
viation obtained from Eq. (|88[) . and relying only on the 
mean string power spectrum, together with the measured 
standard deviation of our bispectrum estimator. There 
is good agreement for both the power law behaviour and 
the overall amplitude. From this section, we conclude 
that although the cosmic string equilateral bispectrum 



is non- vanishing, its measurement at small angles may 
be difficult due to its fast decay compared to the asso- 
ciated standard deviation. This is also the case for the 
squeezed and collapsed configurations, but the formers 
being of higher amplitude, this should be less problem- 
atic. 



IV. COMPARISON WITH DATA 

From both the analytical and numerical results, the 
cosmic string bispectra associated with the flat triangle 
configurations are the best suited to look for string sig- 
natures. However, it is not easy to compare with existing 
bispectrum estimates as much of the literature focuses on 
particular local models of primordial non-Gaussianity in 
which the Newtonian potential is parametrised [42j 



(89) 



The local type of primordial non-Gaussianity produces 
a bispectrum whose maximal amplitude occurs for 
squeezed triangle configurations, as it is the case for the 
cosmic strings. However, as a result of the CMB trans- 
fer functions, a given value of corresponds to os- 
cillating damped patterns of the CMB temperature bis- 
pectrum, which are therefore completely different from 
the power law we have derived for the strings at small 
scales. The current bounds on /^J being precisely ob- 
tained from template matching procedures, they cannot 
be applied to the string bispectrum [It], [28[ . Notice that 
the analysis concerned with the other type of primordial 
non-Gaussianities, such as the equilateral ones charac- 
terised byf^L, are also affected by the CMB transfer 
function [17| and cannot straightforwardly be compared 
to our findings. 

However, a simple calculation we can do at this stage 
is, at a given angular scale, to compare the overall am- 
plitude the CMB bispectrum induced by a primordial 
with the one induced by the presence of cosmic 
strings. For simplicity, let us choose the equilateral con- 
figuration (which is not of maximal amplitude for both 
the strings and local primordial non-Gaussianity). With 
COBE-normalised CMB temperature fluctuations, and 
neglecting Silk damping, Ref. [42j gives the amplitude of 
the oscillating CMB bispectrum induced by primordial 
non-Gaussianity: 



max I B 



(2 x KT 17 )/- 4 !/^! 



(90) 



We can now estimate the size of the equilateral bispec- 
trum from strings if they contribute about 10% of the 
temperature power spectrum at I = 10. In this case 
e ~ 2 x 10~ 5 , and assuming Abelian Higgs strings where 
GU ~ 0.7 x 10~ 6 [3,111, one has 



B, 



-(4xl0- 14 )c 



GU 



£4 A £2 fc 6 ^0.7 X 10 



(91) 
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Matching the two suggests that string fluctuations may 
generate a temperature bispectrum at £ ~ 500 as large 
as the one that would be produced by a local |/nlI of 
around 10 3 ~ 1 is near the peak of the scalar modes 
at £ ~ 500). This number should be considered as an 
order of magnitude estimate since the string bispectrum 
calculation is only accurate for higher £, while Eq. (|5D|) 
is applicable for £ < 1000 [42|. Let us also mention that 
this value for is only valid at I ~ 500 and does not 
represent the value that would be obtained by a local pri- 
mordial bispectrum template matching against the string 
bispectrum. In order to obtain such a value, one should 
perform a Fisher matrix analysis over all the multipolcs 
along the lines of Refs. 0, |45j- Such an analysis is out- 
side the scope of this paper and we leave it for a forth- 
coming work. 



non-Gaussian signals. 

Finally, as the power law dependence of the bispectrum 
variance may suggest, we should expect an even stronger 
non-Gaussian signal from the cosmic string trispectrum 
and we leave its analysis for a forthcoming work. 

Acknowledgments 

This work is partially supported by the Belgian Fed- 
eral Office for Scientific, Technical, and Cultural Af- 
fairs through the Inter-University Attraction Pole Grant 
No. P6/11. The cosmic string simulations have 
been performed thanks to computing support provided 
by the Planck-HFI processing center at the Institut 
d'Astrophysique de Paris. 



CONCLUSION 



APPENDIX A: STRING COORDINATE 
INTEGRATIONS 



We have analytically and numerically derived the CMB 
temperature bispectrum induced by cosmic strings from 
the Gott-Kaiser-Stebbins effect. This effect is the main 
source of string induced CMB anisotropics at small an- 
gle and our result is applicable at those angular scales 
only. The bispectrum is negative for equilateral and all 
isosceles squeezed configurations, having a power law de- 
cay close to £~ 6 . The squeezed isosceles bispectrum is 
significantly amplified for small squeezing angle by a fac- 
tor 9~ 3 . Although the vanishing limit is forbidden by any 
finite detector resolution, one should expect the ongoing 
and future high resolution CMB experiment to be more 
sensitive to these signatures [III [43, [48[ . The collapsed 
isosceles configuration still exhibits the same power law, 
an angle dependency in ip~ 3 (with ip = n — 6) but leads 
to a positive bispectrum. This feature could allow an 
enhancement of the string signal by subtracting the bis- 
pectrum of these two configurations. The numerically 
measured amplitudes have been given in the text and, 
at the particular scale £ = 500, they roughly correspond 
to a I/jjlI — 10 3 f° r the local type of primordial non- 
Gaussianity in the equilateral configurations. This num- 
ber is of illustrative purpose only and should be treated 
with caution. Firstly, the temperature bispectrum in- 
duced by the strings has a different multipole dependency 
than the one coming from primordial non-Gaussianities. 
As a result, the value for /^l quoted above is definitely 
not the one that would be obtain by applying a pri- 
mordial bispectrum template matching over the actual 
string bispectrum. In particular, the current bounds de- 
rived so far on the primordial /nl values do not apply 
for the strings. Moreover, even at the string prediction 
level, there are competing and correlated effects at this 
angular scale, notably the acoustic oscillations that we 
have not considered here. Nevertheless, this value for 
/jJl shows that this is an interesting angular scale to 
search for string non-Gaussianity and that devoted string 
template matching should be used in the search of CMB 



In this Appendix we perform the integrals over the 
string coordinates a a m the expression (|55[) for the bis- 
pectrum. We repeat the expressions for convenience: 

i?(fc 1 ,fc 2 ,fc 3 ) = -6 3 l fcl ^ fc 2 3c f daid<r 2 da 3 

A fcffcpf J (Al) 



x (C ABC D)e-^ D >. 



where 



k lA k 2B k 3c (C ABC D)e-^ = 

J j K 12 [«13ll(£7l3) + K23n(cr 2 3)] V(<T 12 ) 
+K13 [Ki 2 n(cri2) + K23n(cr 3 2)] V(<J3i) 
+ K23 [Kl 2 n(a- 2 l) + Kl 3 n((7 3 l)] V{(J 23 



x exp j-- [ K 13 T(a 13 ) + n 23 T(a 23 ) + K 12 T(a 12 ) } 

(A2) 

Thanks to the symmetry identified earlier, there is essen- 
tially only one integral to do, as the six terms in Eq. |55|) 
in are related by a permutation symmetry. We can define 
a function of three variables 

F(k 12 ,k 23 ,k 31 ) = J dcrido-2dcr 3 n(cri 3 )V(CTi2) 



x exp <{ - [-K 31 r(cr 13 ) - K 23 r(er 2 3) - n 12 T(a 12 ) 



in terms of which 

k A kik% (C ABC D)e-^ D2 > = -{ 



(A3) 



ki 2 k 3 iF(ki 2 ,k 23 , k 3 i) + ki 2 k 23 F(k 12 , k 31 , k 23 ) 
+k 3 iK 12 F(k 31 ,k 23i Ktt) + n 31 n 23 F(n 31l Ki 2 , k 23 ) 
+k 23 ki 2 F(k 23 , k 3 x,ku) + k 23 k 31 F(k 23 ,k 12 , « 3 i)}. 



(A4) 
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A change of coordinates simplifies the integration: 



Hence 



C123 = g(0"l + CT 2 + CT3), 



(A5) 



giving 



■F(kX2>«23)«3x) = y dCTl23dCTi 2 dCT3in(cr 3 i)V(eri2) 
X exp | | [-/«3ir(<73l) - K23r(CT23) ~ KX2F(cti 2 )] } , 



(A6) 



where we have used IT(<73i) = n(o"i3) and T(a 3 i) 
r(o"i3). Inserting the identity 



1 = J da 23 8 (er 2 3 + (T12 + 0-31) 



(A7) 



and using an integral representation for the <5-function 
we have 



31 



F(ki2,k 2 3, «3l) = L J ^ J do-i 2 dcr23dcr. 
x U(a 3 i)V(ai 2 ) exp [iX (cr 23 + a 12 + a 31 )] 



X exp <j -- [Ki 3 r(f73i) + K 2 3r(cr 2 3) + Ki2r(eri 2 )] 



Hence we can factorise 

F(Ki 2 ,K 2 3,K 31 ) = L 

where 

/(A,ki 2 ) = 

J{\ K 23 ) = 



dA 
2- 



(A8) 



Z(A, K12) J(A, K2 3 )^(A, k 3 i), 
(A9) 



do-V F (<r)e iAa -*' ,1 " r(<r) , ( A10 ) 

dcre iAa-iK 23 r(a) ) ( AU ) 



if(A,K 31 ) = J dan^y^-i^^l (A12) 

We evaluate I, J and if for small angles. For n ab large 
and positive, and assuming the relevant limits for the 
correlation functions 

FW^V, V(*)~v 2 , n(a).JV, (A13) 



one gets 



-9 

V Co 

^(ki2, k 23 , «3iJ - 4?ri-p y 



X 



«23 + «12 



(K23K31 + K12K31 + Kl2K 2 3) 3/2 

. . (A17) 

Though we have assumed n a b > in deriving this equa- 
tion, the result is still applicable even if K a (, is negative 
because the combination K23K31 + ki 2 k 3 i + k\ 2 k, 23 in the 
denominator is positive definite. Actually, by using ap- 
proximations (|A13j) . one can derive (|A17|) directly from 
Eq. (|A3|) which is well-defined even for negative n a b be- 
cause of the positive definiteness of the exponent. Defin- 
ing the new function 



.g(Kl2,K23,K 3 l) 



we can write 



K12 + K 23 



(k 23 K 3 1 + K12K31 + Kl 2 K 23 f /2 

(A18) 



B(ki,k 2 ,k 3 ) = -e 3 7rc 



t 4 A ^ k\klkl 



X {«X2«310(k12, «23, «3l) + Kl2K235(Kl2, «3Xj «23) 
+«31«125(«31, «23i «12) + K 3 lK235( K 31, «12, K23) 
+K23«125( K 23, «31i «12) + >«23K3lS( K 23, «12, «3x)}- 

Thanks to g being symmetric in its first two arguments, 
and using k 2 3 + K\ 2 = k 2 (and circular permutations), 
some simplification follows: 



B(ki,k 2 , k 3 ) = -e 3 irc 



v 2 L£ 1 1 

-„2 1,2 7.2 



{ 



i 4 A |2 fcffef fc 3 
/c|K3ig(Ki2, K23, «3l) + fel«23ff(«3X, K12, K23) 
+fcf/«123(« ; 23,'«3X,KX2)}- 

Finally, using the definition of g from Eq. (|A18|) . 



(A20) 



Tf\ \ V / 47r 

i (A, K12) ~ ^ a / exp 

t V K12 



A 2 



«X2* 2 

A 2 



^ 1 / 4?r 

J(A,K23)-T\/ eX P I T 2 

* V «23 V K 23^ 



^(A,k 3 x) 



c 1 



4tt / 2A 2 

1 — I exp 



I K 3 lt 3 V «3X V K 3xi 2 



(A14) 

(A15) 
A 2 



«3xi 2 , 
(A16) 



v 2 L£ 1 1 

B(k 1 ,k 2l k 3 ) = -e 3 irc 



< 4 „4 |2 fc 2 /c 2 fcf 

fc| ^23 + fc|^31 + fc|«12 
(^23K3X + K12K31 + Kl 2 K 2 3) 3/2 



(A21) 
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